Angular invariant quantum mechanics in arbitrary dimension 
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Abstract One dimensional quantum mechanics problems, namely the infinite po- 
tential well, the harmonic oscillator, the free particle, the Dirac delta potential, the 
finite well and the finite barrier are generalized for finite arbitrary dimension in a 
radially symmetric, or angular invariant, manner. This generalization enables the 
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functions, and it also enables connections to multi-dimensional physics theories, like 
string theory. 



INTRODUCTION 



An introductory quantum mechanics course deals with solutions to one-dimensional prob- 



lems, as can be seen in the commonly used textbooks on the subject. Three-dimensional 
problems, like angular momentum, scattering and the hydrogen atom, are not generalizations 
of one-dimensional problems, and they normally require either a specific solution method to 
the Schrodinger equation or an additional symmetry input. As a multidimensional approach 
to one-dimensional problems does not necessarily lead to more relevant multidimensional 
models, it is merely regarded as a curiosity. However, the importance of multi-dimensional 
problems has increased in physics, since string theory has given rise to the possibility that 
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there could be more than three dimensions of space. For example, the Schrodinger equation 
was semi-classically solved in various dimensions in order to quantize pulsating strings jl, 2]. 

On the other hand, the more dimensions a problem has, the more possibilities of motion, 
and the more symmetries it can have to restrict these possibilities. This means that a 
generalization can be a choice, depending of the symmetries of the n— dimensional solution 
of the problem. In this article we answer the question of what happens when the Schrodinger 
equation is solved using a formalism which generalizes one-dimensional problems into angular 
invariant n— dimensional cases. Accordingly, the Schrodinger equation is transformed into 
the Bessel equation, whose solutions appear in physics problems which either have cylindrical 
symmetry or spherical symmetry jsl, 0] . 

The multidimensional approach provides a deeper understanding of the physics of the 
one-dimensional problems, as well as the capabilities and limitations of the mathematical 
apparatus. As some of these results may be scattered throughout literature, it is useful, for 
both students and researchers, to have them presented in a single place. 

This article is organized as follows: in section two, the infinite n— dimensional cylindrical 
quantum well is solved. In the third section, the quantum harmonic oscillator is studied in 
various dimensions, and it is shown that this problem can be described in terms of Bessel 
or Whittaker functions. In section four, the free-particle is analyzed, requiring a particular 
Dirac delta function. Section five deals with the Dirac delta potential. The finite well and 
the finite barrier are dealt with in the sixth seventh sections, respectively. Finally, a brief 
conclusion rounds off this article. 



2. THE INFINITE SQUARE WELL 

The Schrodinger equation can be written schematically as an Eigenvalue equation 

(n 2 + V)V = (1) 

where ^ is the wave-function, £ is the energy, V is the potential and fl 2 is the momentum 
operator. Considering an (n + 1)— dimensional space, the squared momentum operator, 
expressed by means of spherical coordinates, depends on a Laplacian operator that has a 
radial term and an angular term so that 
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where V 2 a=rd is the term of the Laplace operator for the radial coordinate and for the 
angular coordinates. For a radial-only dependent potential, the wave-function, expressed 
as \l/(r, 0) = s(r)w(9), so that w(Q) has n angular variables, splits ([1]) into two equations, 
namely: 

-V 2 r s + e-v-^ = (3) 
s r z 

-Vlw = -M (4) 
w 

where e = 2m£/h 2 , v = 2mV/h 2 and M is a separation constant that is zero for n = 0. 



Equation (j4]) can 
found elsewhere 



3e solved in terms of n— dimensional spherical harmonics, which can be 
6[. In order to solve the radial equation, it has been set 

u{r) 



(5) 



and thus (13]) becomes, 



r V + (n - 2v)r u + [(e - v)r 2 + u(u - n + 1) - M] =0. (6) 

The solution for ([6]) depends on the particular potential v. For an infinite potential well, 
the potential is 

( if r <R 

( 7 ) 

oo if r > R, 

and for a radial dependent-only wave-function, M — 0, the choice v = leads to the 
Bessel equation, 

r 2 u" + ru + (e r 2 — z/ 2 ) u = 0. (8) 
Thus, the wave-function is expressed in terms of Bessel functions 

^«( r ) = ^( aJ A^ r ) +bY u ( v r er )^ (9) 

where a and b are integration constants. The choice of v = 1 ^ would lead to the spherical 
Bessel equation, whose spherical Bessel functions, j u and y u , are related to the usual Bessel 
functions as, 

Ux) = \f^, J v +\ and V»( x ) = \f^ Y v+\' ( 10) 
which generates an identical wave-function, thus, (9) is indeed the most general solution 
to the problem. As Y u (x) is divergent in x = for n > 0, then b = because otherwise 
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the wave-function is not normalizable. The potential is infinite at r > R, thus ^{R) = 0. 
Defining r = as the iV— th zero of J v -i (r), the quantized energy is obtained from 
y/e R = rjy , and it is expressed as 



4 n) = 1 ^ 1 • (ii) 




As R is a free parameter and n is fixed by the geometry, the more excited the level of 
the energy, the more zeros the wave-function in the interval (0, R) has. Normalization is 
obtained from 

J drr^V'er)) 2 = r ^[[j v (V~er)f - J„+i{V~er) J^(V~er)], (12) 
so that J °° dr r n \ \ 2 = 1 implies the normalized wave-function 

*(r) = -J , w , , 2 L (13) 

The wave- function is rotationally invariant. In this sense, the solution calculated above is 
not valid in the n = case; the usual one-dimensional n = solution has anti-symmetrical 
states due to the negative values that the argument of the wave-function has in this case, 
which are not included in a rotationally invariant (n + 1)— dimensional wave-function. 

3. THE HARMONIC OSCILLATOR 

The one-dimensional harmonic oscillator is solved analytically in terms of Hermite poly- 
nomials, as originally demonstrated by Schrodinger [jj. This method uses the asymptotic 
behavior of the wave-function to simplify the problem and to obtain the Hermite equation. 
However, using the variable p = [i r 2 and the index v = in equation (E} it is obtained 



u" + 



2 

1 e (n + l)(n-3) +4M 



u = (14) 



4 4 up 16p 2 

where p = m ^- and the prime means differentiation relative to p. Equation (114p with M = 
is the Whittaker equation, whose general solution is 

u(p) = aM Xiri (p)+nW XiV (p) (15) 
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where A = j- = gfj, rj = ±^^p and a and 6 are integration constants. The Whitaker 
functions M A >r? and Wa,^ m ay be expressed as 

M A ,» = e~V^M(i + ^ - A, 1 + 2r7, p) (16) 

W A) „Go) = e-%p^u(± +77- A, 1 + 277, p) - (17) 

M(p,q,z) and U(p,q,z) are confluent hyper-geometric functions known as Kummer func- 
tions. As always, wave-functions must be normalizable, and Kummer functions diverge if 
the first index p is not a negative integer. Thus, normalizable wave-functions are orthogonal 
polynomials of the order iV e N, obtained through the following relations: 

M(-N, ~ z 2 ) = (-l)«*" H M*) (18) 

M( - N, 2, . 2 ) = (-1)" ^ . ^±1^ (19) 

V ' 2' J K 1 (2iV + 1)! 2z y ' 

V 2 ' 2' / 2^ y ' 
U(-N,a+l,z) = (-1) N (a + l) N M(- N, a + 1, z) = (—1) N N\L^ (z) (21) 



where (a + 1)jv is a Pochhammer symbol, H^(z) are the Hermite polynomials, and L 



(«), 
jv 1 



are the generalized Laguerre polynomials. The one-dimensional solution for the harmonic 
oscillator is obtained by setting n — 0. In this situation, b = in (115j) . according to 
two reasons: when N is even, (12"TT) is divergent in z = and so the wave-function is not 
normalizable, and when N is odd, the wave-function is normalizable but the energy integral 
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£= J drr n \m>\ 2 (22) 

— oo 

is divergent, thus solutions involving (I20I) have to be discarded. The remaining conditions 
(j!8p and ffT9l) imposed on ^ in the n = case, enables us to write the energy spectrum of 
the linear harmonic oscillator 

£ 2N = hw(2N + ^j (23) 
£ 2N+1 = hw(2N+^ (24) 
Using the orthogonality condition for Hermite polynomials 



oo 

J dxe~ x2 H L (x) H K (x) = 2 K K\^8 K , L 



(25) 
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the normalized wave-function for the one-dimensional case is obtained with 



* K = -^=(^Ye-^ 2 H K (^Ir) (26) 



V2 K K\ Vvr, 

where ( )26l) is valid for K = 2N and K = 2N+1. For the n— dimensional case, (12TT) indicates 
that the solution is given in terms of generalized Laguerre polynomials, Lff and thus a = 
is established in (TPol) . without loss of generality. When comparing ( 1T7|) and fT2~Tj) we get 
a = ±^y- • As n > and a > |1, the plus sign must be chosen. Also from (TlTj) and (12T1) . we 
get the energy spectrum 

£^ = ^(2^+— ^— J. (27) 

Using the orthogonality relation 

oo 

J dxx a e- x L$(x)L { «\x) = r(l + a)( N + a y M , N (28) 
o 

the normalized wave-function is 



<Mr) = ^(-1)^ 2r ( iV + 1 \ e Vi^^) (29) 

The result shows that the energy depends explicitly on the angular dimension n and that 
the wave-function has a rotational symmetry, as expected from the angular independence 
imposed by using M = in the Schrodinger equation. 



4. THE FREE PARTICLE 



The equation that describes a free particle is similar to the equation for the infinite well, 
as both have zero potential. The difference resides in the boundary conditions. In the n = 
case, the solution is expressed in terms of complex exponentials, and in the arbitrary n 
case, it is expressed in terms of Hankel functions, which describe cylindrical travelling waves 
expressed in terms of Bessel functions, namely, 

Hl 1 \z) = J u {z)+iY v {z) and H?\z) = J„(z)\iY v (z), (30) (30) 

A general cylindrically symmetrical solution (EJ) of (E]) for the free particle is, 

^ (r) = -4zr (a (V~er) + b H^L (yfer)), (31) 

2 r~ 2~ V 2 2 / 
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in which Hu is a travelling mode towards r = 0, and is a travelling mode towards 
r — 7- oo and a and b are integration constants. As in the n = case, the free-particle 
wave-function is not normalizable and it is understood as a wave packet which obeys, 

oo 

J dr r n ¥ v ^ s = 5 n+1 (77 - 5) (32) 

Q 

where 5 n+1 (i] — 5) is an (n + 1)| dimensional Dirac delta function. By substituting two 
wave-functions like ( 13T1) with energies rj and 5 in ( 1321) . and considering 

H$-\z) = (j_„(z) -e*"J v {z)\ H?\z) = — (- J-vW + e^Mzj),, 

ismuir\ / 2 sin i/7r V / 

and J_ v = (— iyj u , we find that 

00 

5 n+l (e - v ) = V~eJ drrJ v {^~er) J^r) (33) 


which is the definition of a Dirac delta function in terms of Bessel functions, and thus the 
wave-function satisfies the mathematical requirements in order to describe a cylindrically 
symmetric free-particle. However, one physical aspect is missing: the behavior of the wave- 
function at r = 0. There is no external force or internal interaction, thus at this point 
the travelling wave must change direction and maintain intensity. This means that the 
integration constants, which give the wave amplitude of the incoming and outgoing waves, 
must have the same modulus in order to generate equal amplitudes for the wave-function 
at r = 0. As Y v (0) is divergent and cannot contribute to the solution, the wave- function is 
simply, 

WrH^J^Vir). (34) 

2 f 2 2 

Analogous to the n = case, where wave-functions can be expressed in terms of a Fourier 
transform which allows the free particle to be interpreted as a wave packet, the same can 

be achieved here by expressing the free particle wave function as a Hankel transform 

00 

*(r) = J dVe^Ve) j^er) (V~er) 1/2 (35) 


for an appropriate function <f>(y/e). Thus, the analogy between the n = and the arbitrary 
n is complete. Physically the general n— dimensional case is more symmetrical because 
the wave-functions can only have one amplitude for both modes, something that does not 
constrain the n = one-dimensional case. 
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5. THE DELTA FUNCTION POTENTIAL 

In this case there is a Dirac delta function potential, 

V = ±g5(r-R) (36) 

in which g > is the coupling constant of the potential. A negative sign in (136]) means a 
potential well and a positive sign means a potential barrier. Scattered states are possible 
for both signs of the potential, and a bound state occurs in the potential well for negative 
energies, which is discussed in the following subsection. 



5.1. bound state 

The Schrodinger equation with negative energy £ = — \£\ is expressed as 

V 2 * + 7<J(r - R)V = e* (37) 

where 7 = j^g and e = Tpr|£|. The general solution to this problem is given in terms of 
modified Bessel functions I u and K u , integration constants a and b, and v = so that 

tt(r) = ^ (a I u (V~e r) +bK„(V~e r)) . (38) 

The modified Bessel function /„ is divergent at r — >■ 00 and at r — >• 0, ^ — > 00, thus the 
wave-function is, 

f ^ / = a Mv^i) if r<R 

*(r) = < 39 
[^ = 6^ if r >jR . 

At r = i£, ^7 = ^ n and one integration constant is eliminated, 

On the other hand, the first derivative of the wave-function is not continuous at r — R, as 
can be seen from integrating ( 137|) in a r = R neighborhood, which gives, 

A(r n ^') = -7*(i2). (41) 

The e\l/ term is eliminated from (|37|) by integration, and this does not contribute to (). This 
means that the energy sign of the energy is not important in order to determine whether the 
particle is bound or free; all information regarding this is in the potential sign. At r = R, 

A(r n *') = R n (V n (R) - = -R n - v ^(aIu{VeR) - b K„(y/e R)) . (42) 
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Using (??), (jO) . fl42|). and the Wronskian 

K v {x) I u+l {x) + K u (x) I v+1 {x) = -, (43) 

x 

I u (y/eR) K u (^~tR) = — is obtained. (44) 

7 it 

( JHJ) is a transcendental equation and it enables us to determine the energy numerically or 
graphically for each n. However, some particular cases can be calculated. For 

-1 9 9 

x — > oo, KJx) I v (x) — >■ — and thus e = — or 151 = -^tt, (45) 

2x 4 2^ 

which is the unique bound state of this regime. On the other hand if 



1 x 2 

x < 1, then J„(a:) -> — - — — for v>\ (46) 

Iv lv{v A — 1 



so that the energy for this regime is obtained from 



-2^fl-^) (47) 



R \ ~fR 

which is also a one— state solution only, as has been observed in the well-known one- 
dimensional case. 



5.2. scattering state 

In this problem, the particle comes from infinity towards r = and is scattered by a Dirac 
delta well at r = R. In fact, the transmitted wave is totally reflected at r = 0, thus in the 
r < R region the waves have the same intensity in directions; accordingly, the wave-function 
in 

*M =\ / m m \ (48) 

\ = ^(&^ 1} (v^r) + ^ 1) ( v ^r)) if r>i? 

— 1 (2) 

so that ^ = anc i ^ describes the incident wave. From the continuity of the wave 
function and the integration of the Schrodinger equation, we obtain 

a J u - b H { u ] = HP (49) 
a J u+1 + b [-^Hi l) + Hl%) = Hl% - ^Hi 2) . (50) 

All Bessel functions are evaluated at yfeR. Using a Wronskian for Hankel and Bessel 
functions, the above system can be solved for a and b, whose modulus give us the reflection 
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rate and the transmission rate, namely, 

1 C\ 

R=\b\ 2 = l and T=\a\ 2 = 5 ? (51) 

(n 1 RJ u Y u y+ (tt 7j RJ2-2) 2 

A reflection rate equal to one is understandable considering the fact that at r = the wave is 
totally reflected, and as the wave-function describe stationary states, everything coming from 
infinite will be reflected. On the other hand, the transmission rate is something altogether 
more subtle. It may be greater than one, and if it! = 1, it would be expected that T = 0. 
However, there is a reflection of the wave inside the region r < R at r = R, and thus it 
is understandable that, within this region, the intensity of the wave will be greater than 
outside the region; the incoming wave is not immediately reflected to infinity and in fact a 
stationary wave-function is generated by totally reflection at the origin r = and a partial 
transmission at r = R. Thus, T cannot be interpreted as a transmission of the incoming 
wave, but as a relative intensity of the beams in the confined region and open region. Inside 
each region, the incoming beam and the outgoing beam have equal intensities. The intensity 
of the stationary wave drops to zero if the position of the potential R — > 00, but in the zeros 
of the Bessel functions, it is four times greater than the incoming wave, irrespective of how 
far the zero is from the origin of the coordinate system. 

The relative intensity of the wave-functions also enables energy quantization according 
to the value of T, which is an oscillating function. For each particular value of T, there is an 
infinite spectrum of energy where the transmission has this particular value. Thus, it can 
be said that the energy is quantized for this system, because only particular values of the 
energy are permitted. 

One last comment about this case must be made about the wave-function. In the region 
r > R there are incoming and outgoing waves represented by Hankel functions. The intensity 
of these waves is equal, so \b\ 2 = 1, although the coefficients are not necessarily equal, thus, 
by ansatz, b ^ 1. In the r < R region, the situation is different, and the intensity and 
the Hankel functions coefficients functions must be equal in both of the directions because 
Y v — > 00 at r — )• 0, something which does not occur in the r > R region,thus the greater 
generality of the wave-function there. 
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5.3. Dirac delta function barrier 



The case of the Dirac delta function is totally analogous to the scattering of the Dirac 
delta function well tackled above, the only difference being that the sign of 7 sign in the 
potential term is flipped from plus to minus. This change, however, does not alter any of the 
results, which depend only on 7 2 , so the well and the barrier are physically indistinguishable. 

A physical analogy of these models can be executed with a laser beam, which is produced 
from an oscillating light-wave inside a partially reflecting chamber. A light wave is produced 
inside the device and when the intensity of the wave inside it is high enough, the coherent 
light escapes through one of the sides of the chamber, which acts as a barrier, in a situation 
similar to the Delta scattering. To be more realistic, the model would, of course, need a 
source at r = 0. 



6. THE FINITE POTENTIAL WELL 



In this case the potential is 

Vl = if r < R 



(52) 

v. ;c ,• --. /> 

h 2 



vn = -v = -^V if r>R 



where Vo > and the potential describes a cylindrical well whose center is located at 
r = 0. There are two possible solutions: a bounded-state solution with negative energy and 
a scattering-state solution with positive energy. 



6.1. bound states 

This problem has negative energy £ = —\£\ and v — ^ and M = were chosen in (|6]), 
thus obtaining, 

r 2 u" + ru + (Q a r 2 -v 2 )u = Q, (53) 



so that, 



. Qi = v — lei if r < R 
Q a ={ (54) 
Q n = -\e\ if r>R. 
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The general solution to both regions is, 

ui = aJ u (^Qjr) +bY u ( v fQ]r) and (55) 
Ul = cl u (^r) +dK u (^r). (56) 

The wave-function must be finite at r = and at r — > oo, thus b = c = 0. The continuity of 
the wave-function and its first derivative at r — R generates, 



K v {yfeR) J v+1 (yfeR) 









V 




|£| 



(57) 



K v+1 (yfeR) J u {V~eR) 
This transcendental equation is solved numerically for each n, and the intersection points of 
the graphs of both sides give us the quantized energy Quantized energy may be obtained for 
specific cases. If 8 ~ Vq, then J v — > 0, thus the quantized energy comes from, \/~QiR = 



= (fl. («) 



2 



where is the iV|th zero of J v . In this situation R ^> x^ . If n = 0, Bessel functions 
turn to trigonometric functions and the known regular spacing among these zeros for the 
one dimensional well appears. Other quantizing possibilities come from the Vo S> £ regime. 
In this case the quantized energy comes from the zero of J u +i, namely y/eR = x^ . This 
case can be understood as a deep well or tiny energy. 



6.2. scattering states 

In this case the state has positive energy and the potential is the same used in the bounded 
states (152]) . As discussed in the case of the free particle, the intensity of the wave-function 
is maintained at r = 0, and so the solution is simply, 



^ I = a h — - H h — L if r<R 



and 



*(0 = { r (Z, (59) 

*ii = b J -^± if r<R 

so that V = e + vq. From the continuity of the wave-function and its first derivative, we 
obtain \b\ 2 = 1 and, 

rp | ,2 16/(7T6i? 2 ) 

( Jv Jv+l — A/ -£ + 1 Jv Jv+1 J + ( Jv Yv + l — \l -£ + 1 Yv Jv+1 
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so that J v = J v {yV R), J u = J u (y/eR) and Y v = Y v (^feR). The result is compatible 
with the situation for scattered states observed in the delta function model, where the 
incoming wave is totally reflected at r = R because it is totally reflected at r = 0, and 
the intensity of the wave-function is greater in the r < R region because the outgoing wave 
is partially transmitted at r = R. For each particular value of T, which is an oscillating 
function, there is an infinite energy spectrum that gives us this value, and then the energy 
is quantized according to this particular value. This result has also already been obtained 
for the scattering in the Dirac potential. 

7. CONCLUSION 

In this article several one- dimensional quantum mechanics problems have been generalized 
in an angular invariant manner. The results confirm expectations such as a contribution to 
the zero-point energy of the harmonic oscillator due to the angular dimension, and several 
results are not so obvious: the equality of the wave-function intensities to free particles, the 
existence of various quantized states in the Dirac delta potential and the quantum scattering 
states for the finite well. It is hoped that this set of results will be useful in understanding 
quantum-mechanics problems that link angular invariance and multi-dimensionality. 
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